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Abstract 

We study the addition of quantum fermionic impurities to the A/" = 6 supersymmetric 
Chern-Simons-matter theories in 2 -|- 1 spacetime dimensions. The impurities are introduced 
by means of Wilson loops in the antisymmetric representation of the gauge group. In a 
holographic setup, the system is represented by considering D6-branes probing the AdS^^ x 
CP'^ background of type IIA supergravity. We study the thermodynamic properties of the 
system and show how a Kondo lattice model with holographic dimers can be constructed. 
By computing the Kaluza-Klein fluctuation modes of the probe brane we determine the 
complete spectrum of dimensions of the impurity operators. A very rich structure is found, 
depending both on the Kaluza-Klein quantum numbers and on the filling fraction of the 
impurities. 
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1 Introduction 



The gauge/gravity correspondence [H [2] has been an invaluable analytic tool which has 
provided new insights on the the dynamics of gauge theories in the strong coupling regime. 
Indeed, in the framework of string theory, the gauge symmetry is realized as the worldvolume 
symmetry of a stack of coincident branes and the corresponding gravitational background 
encodes relevant properties of the gauge theory in the 't Hooft limit. The richness of string 
theory allows several possibilities to deform the theory by adding extra degrees of freedom. 
In particular, one can add new sets of branes which intersect with the color branes that host 
the gauge symmetry and create a defect in the gauge theory. By studying the dynamics of 
the extra branes in the original background one can extract useful dynamical information of 
the corresponding gauge theory with impurities. 

In this paper we will concentrate on analyzing the case in which the defects are point- 
like. The canonical example of this type of systems is A/" = 4 Super- Yang-Mills (SYM) with 
fermionic impurities. The brane realization of this system is obtained by adding D5-branes 
to the AdS^ x background, dual to A/" = 4 SYM. The D5-branes are extended in the 
holographic direction and wrap an S'^ inside the in such a way that their worldvolume is 
an AdS2 x §^ submanifold of AdS^ x [31 H] . The that the D5-branes wrap are obtained by 
fixing the latitude inside the §^ to a discrete set of angles, which are obtained by imposing 
a quantization condition These configurations can be regarded as describing bound 
states of fundamental strings extended along the holographic coordinate. Equivalently, they 
represent Wilson lines in the antisymmetric representation of the SU{N) gauge group of the 
N' = 4 SYM theory [5]. As shown in detail in [6], these Wilson loops can be obtained by 
integrating out a fermionic impurity, which is localized in the spatial gauge theory directions 
and whose occupation number is just the number of strings carried by the D5-brane. More 
recently, these configurations have been utilized [71 [8] to construct holographic dimer models, 
by considering D5-branes which connect two spatially separated impurities on the boundary 
of AdS^. It has been also argued that this setup of D3- and D5-branes realizes holographically 
the maximally supersymmetric Kondo model [9l [10] (see also [TTl [12]). 

The holographic study of gauge theories with impurities may be useful to understand 
the applications of the AdS/CFT correspondence to condensed matter theories. For this 
reason it is very important to have examples of such systems at our disposal, in which both 
sides of the correspondence are clearly identified and one can perform non-trivial tests of 
the duality. Following this line of thought, we study in this paper the addition of fermionic 
impurities to the Af = 6 Chern-Simons-matter theory proposed by Aharony et al. (ABJM) 
in ref. [13] (based on previous results by Bagger and Lambert [H] and Gustavsson [T5]). 
The ABJM theory describes the dynamics of multiple M2-branes at a C^fZk singularity 
and is a U{N) x U{N) (2+l)-dimensional Chern-Simons gauge theory with levels {k,—k) 
and bifundamental matter fields. In the large limit this theory admits a supergravity 
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description in M-theory in terms of the AdS4^ x S''/Zfc geometry with flux. Moreover, when the 
Chern-Simons level k is large the system is better described in terms of type IIA supergravity. 
In this ten-dimensional description the geometry is of the form 74^5*4 x CP^ with fluxes of 
the RR two- and four- forms. 

It was proposed in ref. [16j that the Wilson lines in the antisymmetric representations 
of the gauge group of the ABJM theory are holographically realized as D6-branes extended 
along an AdS2 C AdS^ and wrapping a five- dimensional submanifold of CP^. These are, 
therefore, the natural brane configurations to consider in order to construct the holographic 
dual of the Chern- Simons-matter theory with fermionic impurities. In this paper we will 
confirm this identification and we will study in detail some of the properties of the sys- 
tem with impurities. This brane setup constitutes the holographic description of a Kondo 
model in which the impurities are coupled to an ambient interacting CFT in three spacetime 
dimensions. 

In the configurations we will focus on, the probe D6-branes capture part of the RR flux 
of the background and, as a consequence, an electric worldvolume field must be turned on. 
This worldvolume field can be eliminated from the lagrangian of the probe by means of a 
Legendre transformation, after using the quantization condition found in [1]. The result 
is an effective one-dimensional problem which can be solved by means of a first-order BPS 
differential equation which depends on the quantization integer n (that represents the number 
of fundamental strings). In this paper we will mostly center in studying the simplest solutions 
of the BPS equation, which give rise to flux tubes made of n strings with < n < (A^ is 
the rank of the U{N) gauge groups). There are, however, more general solutions which are 
related to the baryon vertex of the ABJM theory. A first analysis of these baryonic solutions 
is relegated to an appendix. We will leave a more complete study for future work. 

In the flux tube solutions used to create the impurities, the embedding of the D6-branes 
is such that one of the angular coordinates that parametrize the CP'^ is fixed to a constant 
which depends on the number n of strings. In the dual field theory n is just the number of 
defect fermions and the ratio n/N is the filling fraction of the impurities. The D6-branes 
wrap a squashed T^'^ space inside the internal CP^. Recall that the T^'^ space is a five- 
dimensional Sasaki-Einstein space whose metric can be represented as a U{1) fibration of 
X S^. In our case this space is internally deformed by squashing factors which preserve 
the fiber structure and depend non-trivially on the filling fraction of the impurities. 

Once the setup required to add impurities is clearly identified, one can study many prop- 
erties of the system. We will first analyze its thermodynamic properties by considering a 
probe brane in a black hole background dual to the ABJM theory at finite temperature. 
We will also study the configurations in which two flux tubes ending on two different points 
on the boundary of AdS^ are connected in the bulk and we will show that there exists a 
dimerization transition similar to the one described in refs. [71IB|. In order to establish the 
stability of the brane configuration, we will analyze in detail the fluctuations of a brane probe 
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at zero temperature. We will be able to decouple the different modes and read the conformal 
dimensions of the various dual operators. To carry out this analysis we will have to develop 
the harmonic calculus on the squashed T^'^ space. As a result we will uncover a rich structure 
which depends both on the filling fraction of the impurities and on the quantum numbers of 
the Kaluza-Klein harmonics of the squashed T^'^. As happens in other conifold-like theories, 
the conformal dimensions will not be, in general, rational numbers [17]. We will be able to 
characterize the normalizable modes which do correspond to impurity operators of rational 
dimensions. These operators should be the ones protected by the supersymmetry of the 
system. 

The rest of this paper is organized as follows. In section [2] we briefly review the gravity 
dual of the ABJM theory in the type IIA supergravity. In section [3] we study a class of 
configurations of the probe D6-branes in which they wrap a co-dimension one submanifold 
of CP^. Section m is devoted to analyzing the flux tube embeddings. The thermodynamic 
properties of these configurations are studied in section [5] by considering embeddings in the 
black hole version of the ABJM metric. In this section we compute holographically the free 
energy and entropy of the impurity. In section [6] we study D6-brane embeddings in which 
the brane is hanging from the boundary and ends on two different points. In this section we 
compute the free energy of this hanging tube configuration and we explore the corresponding 
dimerization transition. In sections [7] and |8] we analyze in detail the fluctuation modes of 
the D6-branes around the flux tube and dimer configurations at zero temperature. We check 
their stability and compute the conformal dimensions of the dual impurity fields. In section 
[9] we summarize our results, extract some conclusions and point out some lines for possible 
future research. The paper is completed with several appendices. In appendix |A] we briefly 
present the BPS worldvolume equation and find its solutions. In appendix [B] we derive the 
lagrangian for the fluctuations of the probe brane. In appendix [C] we study the laplacian 
for the squashed T^'^. In appendix [D] we develop the harmonic calculus for these spaces. 
Finally, in appendix |E] we review some basic facts about Heun's equation and its connection 
to a quantum integrable model. 



2 The ABJM background 

The ABJM background [13] of type IIA string theory is characterized by a ten-dimensional 
metric of the form AdS^ x CP^: 

ds' = -^{ds\,s, + 4rf4p3) , (2.1) 

where ds\^g^ and ds^-^-j. are, respectively, the 74^5*4 and CP^ metrics. The former can be 
written, in Poincare coordinates, as: 

ds\,s,=r^{-de + dx- + dy^) + '^4- (2-2) 
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In (12. ip we have represented the AdS^ radius in terms of R and k. Actually, the ABJM 
solution depends on two integers N and k which correspond, in the dual gauge theory, to 
the rank of the gauge groups and to the Chern-Simons levels, respectively. The AdS^ radius 
in (12. ip can be written in terms of these two integers as: 



R^ 2N , , 

where we are taking the Regge slope a' = 1. The metric cis^pa in (12.10 is the canonical 
Fubini-Study metric. In order to write its expression in a convenient form, let us consider 
four complex coordinates [i = I,-- - ,4) such that Yli — 1- Therefore, these z^^s 
parametrize a seven-sphere. We will represent them as: 

1 a 01 i(2ipi+x + i) 2 

z = cos — cos — ^ ^^'^ , z 
2 2 

z' = sin- sin^e^(2^^-^ + «) , z' = sin- cos ^ ("^^^ -^ + «) . (2.4) 
2 2 ' 2 2 ^ ^ 

It is straightforward to verify that the metric of the seven sphere can be written in these 
coordinates as a U{1) bundle over CP^, namely: 

dslr = dslp3 + {dC + Ay . (2.5) 

In our parametrization C ~ ^ the CP^ metric (is^ps, as well as the U{1) connection A, 
are given by: 



a 


. e. 


cos- 


sm- 






. a 




sm — 


COS- 


2 





ds' - ^ 



da'^ + cos^ ^{d0l + sin^ 0i d(fl ) + sin^ ^ {d0l + sin^ 02 difl ) + 

.nO; 2*^/7 /IT /ir\5 

+ sm — cos — [dx + cos0id(pi + cosc^2"V^2j 

A = cosadx + 2 cos^ ^ cos0idipi — 2 sin^ ^ cos02d(p2 . (2.6) 

The ranges of the angles are < a,0i,02, < vr, < fi,'^2 < 27r and < x < 47r. Notice 
that, in the representation (12. 6p . the submanifolds with a = constant are squashed T^'^ 
five- dimensional spaces, which foliate the entire CP^. 

The ABJM type IIA background is also endowed with RR two- and four- forms F2 and 
F4, as well as a constant dilaton 0, which are given by: 

6^'^ = ^, F, = ^R'nAds,, F2 = ^dA, (2.7) 

where flAdSi is the volume element of the AdS^ metric (12. 2p . This background is a good 
supergravity dual of the Chern-Simons matter theory when the AdS4^ radius is large in string 
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units (and the corresponding curvature is small) and when the string coupling e*^ is small. 
From (12.31) and (\2.7\\ it follows [13] that these conditions are satisfied when << k << N. 

In type IIA supergravity the RR six-form Fq is defined as the Hodge dual of the RR 
four- form F4. For the ABJM background Fg is given by: 

^6 = * = ^ siir' ada A €5 , (2.8) 

where €5 is the five-form: 

€5 = sin 9i sin 62 dOi A A dx A dipi A d(p2 ■ (2.9) 

Let us now find the five-form potential C5 for Fq {i.e. such that Fq = dC^). This potential 
will be very useful in the analysis of the probe D6-branes of section [31 First of all, we define 
the function C{a) as the solution of the following first-order equation: 

^ = -3(sina)' , C(0) = . (2.10) 

da 

This function is just: 

C(a) = cosa(sin2a + 2) - 2 . (2.11) 
Then, it is immediate to verify that the RR potential five-form C5 can be taken as: 

C, = -^C{a)e,. (2.12) 



3 D6-brane probes 

Let us consider a D6-brane extended along the radial direction r and wrapping a five- 
dimensional submanifold inside the CP^ at fixed values of the Minkowski coordinates x,y. 
Notice that a configuration of this type creates a point-like defect in the gauge theory. In 
order to describe it in detail, we will take the following worldvolume coordinates: 

C = {t,r,^') = {t,r,ei,^i,e2,f2,x) ■ (3.1) 

We will consider embeddings in which the angle a is a function of the radial coordinate r, 
namely: 

a = a{r) . (3.2) 

Notice that the hypersurface a = constant defines a co-dimension one submanifold of the 
internal CP^. In addition, we will switch on an electric worldvolume gauge field Fq^. Actually, 
such a worldvolume field is sourced by the RR potential C5 through the Wess-Zumino (WZ) 
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term of the D6-brane action (see below). The Dirac-Born-Infeld (DBI) piece of the brane 
action is just: 

Sdbi = -Td, j d'Ce-* V-det(^7 + F) , (3.3) 

where g is the induced metric on the worldvolume of the D6-brane and Tg is its tension. The 
DBI determinant, integrated over the five angles 7* in ( 13. ip is: 



rf^7e-V-det(^ + F) = -^^-^^ (sine)' ^1 + («f - Fl , (3.4) 

where d^'y = dOi A d62 A dx /\ dipi A dip 2 and a' denotes da /dr. Let us define the effective 
DBI lagrangian density Cdbi as: 



Sdbi = / dtdr Cdbi ■ (3.5) 

Then, one has: 



^DBi = - ^lo^sp (sina)' ^J 1 + r^af - ^ For ■ (3-6) 
Moreover, the WZ action is: 

Swz = n Jc^AF = J drdtCwz • (3.7) 

By using the expression of C5 written in f l2.12p . one arrives at the following expression of 
Cwz'- 

Cwz = -^^C{a)For. (3.8) 
Then, the total lagrangian density C = Cdbi + Cwz can be written as: 



C 



^ (sina)^ ^1 + r2 {af - ^ Fl + C{a) For 



(3.9) 



The equation of motion of the gauge field derived from the lagrangian density (13.91) implies 
that: 

dC 

constant , (3.10) 



dFor 

which is nothing but the Gauss law. By using the quantization condition derived in |4] , one 
can relate this constant to the number n of strings (quarks) of the flux tube, namely: 

H ^ "T, . (3.11) 
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where Tj is the tension of the fundamental string and n G Z is the fundamental string charge 
carried by the D6-brane. For our lagrangian density (13.91) we get: 



dC 



dFf 



Or 



28 TT^k 



4k 



For ( sin a)' 



l + r2 {a'Y - 



C{a) 



Or 



(3.12) 



By using this expression in (I3.1ip we can determine F^r in terms of the other variables. 
Indeed, let us define the function Cn{a) as: 



4r? 

C„(«) ^ C{a) + - , 



whose explicit expression is: 



Cn{ce) = cosa(sin^a + 2) — 2 



4n 
iV 



Notice that Cn{a) satisfies the equation: 

dCJa) 



da 



3 ( sin a)' 



4ri 



(3.13) 



(3.14) 



(3.15) 



By combining (I3.1ip and (I3.12p one can get a expression of For in terms of the function 
Cn{a), namely: 



Or 



i?3 Jl + r^a'Y 



Ak 



( sin a) + Cn{ci) 



CJa) 



(3.16) 



In order to eliminate the electric field from the equations of motion, let us next compute 
the hamiltonian of the system by performing the Legendre transformation of the lagrangian 
density ([33]): 



H 



dr 



dC 



Or 



dF, 



C 



Or 



drH . 



By explicit calculation, one gets: 



Or 



dFf 



- £ 



1 + r2 {af 



Or 



210 ^3^2 



( sin ay 



(3.17) 



(3.18) 



Moreover, from the expression of Fq,. written in (I3.16p . we deduce: 

(sina)^ \/ {smaf + CJ 



a 



16fc2 77-2 
ij6 ^Or 



1 + r2 [a'Y 



(3.19) 
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Thus, the hamiltonian density "H can be written as: 

^ = ^TS^P ^1 + ^' \/(sin«)' + C„(a)2 . (3.20) 

The different embeddings of the D6-brane described by our ansatz can be found by integrating 
the Euler-Lagrange equation derived from fl3.20p . In the main part of this paper we will 
concentrate on studying the simplest of these configurations, namely those for which a is a 
constant. These embeddings will be analyzed in detail in following sections. In appendix [XI 
we will present embeddings with non-constant a, which satisfy a first-order BPS equation. 
These last configurations are related to the baryon vertex of the ABJM theory (see |1H1 [19] 
for a similar analysis in the AdS^ x background). We will leave their detailed study for 
a future work. 



4 Flux tube configurations 

Let us now look for solutions of the equations of motion which have a constant a angle. 
Clearly, we should require that: 

^« " (4.1) 



da 



. 



a'=0 



Therefore, we need to compute the derivative of the second square root in (I3.20p . With this 
purpose in mind, we define the following function: 

, 2 



As: 



we get: 



d_ 

da 



m 

da 



A„(a) = (sin a) cos a — C„(a) . 

= 6 ( sin a)'^ A„(a) 
( sin a)^ 



( sin a)^ + Cn{a 



_ 3R^ 
-0 ~ 2i0 7r3A;2 



AJa) . 



(sina)*' + Cni 



(4.2) 
(4.3) 

(4.4) 



a 



Therefore, the non-trivial configurations with constant a occur when a = an, where is a 
solution of the equation: 

A„(a„) = . (4.5) 

Let us now study the solutions of this equation. By plugging the definition of C„(a) written 
in fl3.14p into (14.21) . one can obtain the explicit expression of A„(a), which is given by: 



AJa) 



N 



cos a 



Thus, the roots of (14.51) are 



cos ttr 



1 - 



2n 
iV 



(4.6) 
(4.7) 
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Notice that there is a unique solution a;„ G [0, vr] for every n in the range < n < N . This 
configuration was proposed in reference [16] as the holographic dual of the Wilson loop of 
the ABJM theory in the antisymmetric representations of the gauge group. The fact that 
n < N supports this identification and is a manifestation of the so-called stringy exclusion 
principle. The energy density for these configurations is just: 

jY ^3 2 

En = 'H{a = an) = ^ -rr (sinaj . (4.8) 

Using the value of written in (14. 7p . we get the following binomial law for the energy 
density: 

where we have written the result in terms of the 't Hooft coupling A of the Chern-Simons- 
matter theory (A = N/k). Notice that En is invariant under the change n — )■ — n, as it 
should for the holographic dual of an object transforming in the anti-symmetric representa- 
tion of the gauge group. Notice also that for small n (or large N) the energy density En is 
just equal to n times the tension of a fundamental string extended along the radial direction 
in the geometry (12. ip (which is given hjTfR^/Ak = ^f\j2). This is, again, in favor of the 
interpretation of this configuration as a bound state of fundamental strings, as those studied 
in |1] for backgrounds generated by Dp-branes. The fact that En < nTf E? /Ak is indicating 
that the formation of the bound state is energetically favored and that the configuration is 
stable. We will check explicitly this stabihty by computing the fiuctuation spectra of the 
brane in section [71 It is also interesting to point out here that, as checked explicitly in [16] . 
these configurations are 1/6 BPS, i.e. they preserve four real supercharges. 

The constant electric field /or corresponding to this configuration can be obtained from 
(I3.16P by substituting a by a„, namely: 

for = cos an ■ (4.10) 
Notice that the induced metric on the D6-brane worldvolume takes the form: 



3 



4k 



(4.11) 



which is of the form AdS2 x T ' , with T ' being a squashed T ' space, with metric: 

d4,,i = gijd-f'd-f' = cos^^{del + sin" 01 difil) + sm^'^{del + sin'^a^^V^a) + 



-t- sin^ cos^ ( (ix + cos 6'i(iv9i + cos6'2(iv^2 ) • (4.12) 
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Taking into account that the angles a„ in (14. 7p satisfy: 



sm - = - , ^os - = , (4.13) 



it follows that the metric ds"^-^ ^ can be written in the form: 



ds|i,i = ^ ( d9l + sm^ c/(^2 j ^ _ ^^2 ^ gj^2 ^^2 j ^ 

+ — - ( rfx + COS 9id(pi + cos ^2 rf<^2 ) • (4.14) 

Notice that the line element in fl4.14p collapses to the one of to a two-sphere when n = 0,N. 
This implies that the corresponding D6-brane configurations are singular and, accordingly, 
we will assume from now on that < n < N. Moreover, the metric f l4.14p displays the 
following symmetry: 

(^1,(^1) ^ (^2,^2), n ^ N-n. (4.15) 

The configuration of D6-branes described above creates a point-like defect in the bulk Chern- 
Simons matter theory. The n fundamental strings carried out by the D6-branes introduce n 
fermions localized on the defect, which can be regarded as fermionic impurities in the Chern- 
Simons-matter theory. As pointed out in [6] for the similar configurations in A/" = 4 SYM, 
by integrating out the defect fermions one should get an effective theory for the bulk fields in 
which there is an insertion of a Wilson loop in the antisymmetric representation of the gauge 
group. This argument establishes a correspondence between the Wilson loop interpretation 
of these configurations (as proposed in [16]) and the model of quantum impurities developed 
in this paper. Notice that our setup describes a Kondo model in which the ambient theory 
is an interacting supersymmetric Chern-Simons-matter theory. The ratio: 

^ = ^ , (4.16) 

will be referred to as the filling fraction. Notice that the metric fl4.14p depends on u and 
that under the symmetry fl4.15p the filling fraction transforms as z/ — 1 — z/ and, therefore, 
f l4.15p can be regarded as the realization, in our holographic model, of the particle-hole 
transformation. It is also clear from ( I4.15P that this particle-hole transformation is realized 
geometrically by the exchange of the two S^'s inside the T^'^. The properties of the defect 
operators depend on u and, in our setup, they will be encoded in the squashing of the internal 
j'l,! geometry. 

In what follows we will study holographically the properties of the impurity model. We 
will start in the next section by analyzing its basic thermodynamical properties. 
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5 Impurities at finite temperature 



Let us now consider the ABJM theory at finite temperature. The corresponding background 
is obtained from the supersymmetric one by substituting the AdS^ factor by the geometry 
of a black hole in AdS4, namely: 



with: 



ds^ = -^(^sIh^ + Adsl^s) , (5.1) 

dsln, = -r' f{r) dt' + + [ dx' + dy'], (5.2) 

where /(r) is the function: 

fir) = 1 - ^ , (5.3) 
and vq is related to the temperature of the black hole by means of the expression: 



27r 



1 d 



^ dr 



v^)] (5.4) 

/ J r=r 47r 



In order to compute the free energy and entropy of the fiux tube configuration described 
in section HI let us evaluate the Euclidean action of one of such configurations that extends 
from the horizon at r = tq until a cutoff value of the radial coordinate r = r cutoff- The result 
is: 

En 

Id6 = (rcutoff-ro) , (5.5) 

where En is the tension written in (14. 9 p and we have integrated over a periodic Euclidean time 
circle of period 1/T. The action (15. 5 p is divergent as r cutoff and must be renormalized 

by subtracting a counterterm. Following the holographic renormalization approach [201 
[2T] . the counterterm is just the action of the brane embedded in the metric (15. ip at zero 
temperature and extended from the origin at r = until r = Vcutoffi namely: 

E 

T counter _ r ('^fi^ 

^D6 - JTJ ^cutoff , 10. bj 

where T' is given by: 

J_ _ fjf^ cutoff) j-g 

and /(r) is the function defined in (15. 3p . To evaluate (15. 6p we have chosen the time period 
1/T' in such a way that it corresponds to the same induced metric at r = r cutoff as the 
one of the brane embedded in the black hole metric (15.11) [7j. Therefore, the renormalized 
euclidean action is just: 

rrenor mail zed i- / t TCOunter\ En q\ 
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In terms of this renormalized action, the impurity free energy is defined as: 

pimp rp jrenormalized 

Therefore, by using the value of the renormahzed action, we get: 

F'^P = -En ro . (5.10) 

If one ignores non-abelian interactions, the free energy for M D6-branes is just M times the 
result (15.1 op . By using (14. 9p and (15. 4p . one gets in this case: 

pimp ^ _W2 yXnMf 1 - . (5.11) 



(5.12) 



3 \ N 

The impurity entropy is related to the free energy as: 

gpimp 



-•imp 



dT ■ 

In our case, we get: 

^''"P = v^nM(l - z/) . (5.13) 

3 

Notice that the impurity entropy (I5.13P depends on (square root) of the 't Hooft coupling A, 
as also happens for the Af = 4 impurities [3 HO] . This fact is a consequence of the interacting 
nature of the ambient CFT. As a function of the filling fraction u = n/N, the entropy (15.130 
is maximized at the half-filling point u = 1/2, which is the fixed point of the particle-hole 
symmetry. This behavior is similar to the one found in the N' = 4 theory [TJ |10]. Notice, 
however, that the dependence on u in our case is analytic. 



6 Hanging flux tubes 

In this section we will consider configurations in which the flux tube does not reach the origin 
of the holographic coordinate but instead it starts from the boundary of AdS^, reaches a 
minimal value of r and comes back to the boundary. We will see that, in order to have 
such a hanging string configuration, one must have a non-constant cartesian coordinate 
for the embedding. Therefore, these connected configurations give rise to dimers on the 
boundary field theory, i.e. to flux tubes connecting n quarks and n antiquarks separated a 
finite distance. 

Let us therefore consider a D6-brane probe as in section [3] and let us choose the same set 
of worldvolume coordinates as in (13.11) . We embed the D6-branes in the black hole geometry 
(15. ip and we consider an ansatz in which both a and the cartesian coordinate x vary with 
r, namely: 

a = a{r) , x = x{r) . (6.1) 
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The DBI lagrangian density is now given by: 



^DBi = -^T^^ (sm«)^ ^l + rV(r) {af + r^W (x')^ - . (6.2) 

The WZ lagrangian is still given by f l3.8p . Thus, the total lagrangian density is: 



C 
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— (sin a)' ^1+ r2/(r) (af + r^W (a;') - -^^ F^. + C (a) For 



(6.3) 

We can now eliminate For following the same steps as in section |3] [i.e. by integrating the 
Gauss law and by performing a Legendre transformation). Clearly, the result that is obtained 
in this way is the same as in section |3l after performing the substitution: 

r^{af r^f{r) («')' + f{r) {x'f . (6.4) 

Therefore, the worldvolume electric gauge field is given by: 

i?3 V^l + rV(r) (gf + rV(r) {xf 

For = ^ / g = Cn{a) , (6.5) 

J ( sin a) + Cn{ay 

where Cn(a) is the same function as in (13.141) . After the Legendre transform we now arrive 
at the following routhian: 

^ = 2^^^ \ll + r^f{r) {a'f + r'f{r) {x'f ^( sin a)' + C.(a)2 . (6.6) 

It follows that in this case there still exist solutions with constant a for the same angles 
a = a„ as in (14. 7p . From now on we will restrict ourselves to the case in which a = In 
this case: 

U = v/l + r4/(r) {x'Y , (6.7) 

where E^ is just the energy written in (14. 9 p for the a = a„ configuration with x constant. 
The equation of motion of x{r) has a first integral given by: 

where A is a constant of integration. Then, one has: 

K , ^'^^^^^^ = A . (6.9) 
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Clearly, having A = in (16. 9 P means that x{r) is not constant. From this expression we can 
obtain x' as a function of r: 

A 1 , , 

x' = ± ^= , (6.10) 

^V/(0 ^r4/(r) - A2 



where A is defined as: 



Notice that, as: 



A ^ A . (6.11) 



^l + rV(r) {xf = /'^^ , (6.12) 

^r4/(r) - A2 

the electric worldvolume field Fq^ is given by the following function of the radial coordinate: 



i^or = yr , = cos an . (6.13) 

V/W - A2 



It is important to point out that the two signs in (I6.10p are needed in order to describe 
a hanging D6-brane that comes from r = oo until a minimal r and then comes back to the 
boundary. The actual minimal value of r, i.e. the radial coordinate of the turning point, is 
the solution r* of the following quadratic equation: 

- r^r, = A^ . (6.14) 

By integrating the equation fl6.10p for x{r), we get: 

x{r) = ±A [ . (6.15) 



Let us introduce in this integral a new variable z, defined as 2; = p/tq. One gets: 

= ±1 (6 16) 

where q is defined as: 

g ^ ^ , (6.17) 
'0 

and z^: is the root of the following quartic equation: 

z, -q^ = . (6.18) 
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Notice that z^, = r^/r^ and that q parametrizes, through the relation f l6.18p . the radial 
position of the turning point. Moreover, when g = eq. fl6.18p can be solved to give = 1 
(or = To). Therefore, it follows from fl6.17p and fl6.10p that x is constant for q = 0, which 
implies that these configurations correspond to branes that go straight from the boundary 
to the horizon of the black hole. Coming back to the general connected configuration, let us 
suppose that the brane starts at x = —L/2 at r = oo and comes back to the boundary at 
X = +L/2 {i.e. L is the separation of the n-quark and the n-antiquark) . Then from (16.160 
one gets that L as a function of q is given by : 



2q 



dz 



^{z^ -z){z^-z- g2) ■ 



(6.19) 



In figured] we have plotted r^L = ^ T L as a function of q. It is clear from this figure that 



Separation of the n-quark and the n-antiquark 
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Figure 1: In this figure we plot Ltq as a function of q, as given by (I6.19p . The maximal 
value of Ltq occurs when q ^ 0.779. For this value of q we get LtqK. 0.906. 

there exists a maximal value of L for every temperature T and that for fixed L there are no 
connected configurations for large enough temperature. Moreover, when ro L is small enough 
there are two possible values of q which, according to (I6.18p . correspond to two possible values 
of r^,. In order to determine which of these two configurations is the dominant one, let us 
compute the free energy as a function of q. First of all, notice that plugging eq. (I6.12p into 
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the expression of the energy density given by (16.71) . we get: 



(6.20) 



Let us integrate this energy density for a connected configuration whose ends are located at 
r = rmax- One gets: 



E = 2En 

^Jp'f{p)-A' 
Changing variables in (16.211) from p to z = p/ro, we get: 



dp . 



(6.21) 



E = 2Er.ro 



z* - z 



dz 



(6.22) 



This integral is divergent when Zmax oo. In order to regulate it and get the free energy, 
we follow the procedure in [221 123] and subtract the energy corresponding to a disconnected 
configuration that reaches the point r = Tmax- This energy is: 



2E„ 



dp = 2Enro 



dz + z^ 



(6.23) 



The regulated free energy for the connected configuration (after sending z^a.^ oo) is given 
by: 



2E„ro 



1 ] dz — z^ 



(6.24) 



z^ — z — q'^ 

Notice that for g = (and 2;* = 1) the free energy in fl6.24p F = — 2i?„ro, i.e. twice the 
result in (IS.lOp and (15. lip , which confirms that this q = configuration corresponds to two 
straight disconnected branes stretching from the boundary to the horizon. The free energy 
F as a function of q has been plotted in figure [2l By comparing figures [1] and [2] it is easy 
to verify that, when tq L is small enough to allow two solutions for q, the configuration with 
higher q [i.e. the one that penetrates less deeply in the bulk) is the one that has lower free 
energy and, therefore, is the dominant one. It is also clear from figure |2] that, for low enough 
temperature, this configuration has a free energy lower than —2En and, therefore, it is 
also dominant over the disconnected one. This reconnection of branes at low temperature is 
the dimerization transition described in [7]. Its characteristics are very similar to the ones 
of the TV = 4 four-dimensional theory and will not be elaborated further here. 



7 Impurity fluctuations 

Let us now consider fluctuations around the static configurations studied above. We will 
allow to fluctuate the coordinates transverse to the D6-brane (namely the angle a and the 
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Free energy for hanging flux tubes 
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q 

Figure 2: Plot of F/ (2£'„ro) as a function of the parameter q. The curve reaches a maximum 
at g ~ 0.789, where it takes the value F/{2Enro) ~ —0.923. The dimerization transition 
starts at g ~ 1.747, where F = — 2i?„ro. 



cartesian coordinates x and y), as well as the worldvolume field F. Thus, we write: 

«n + ^, F = f + f , X = X + X , (7.1) 



a 



where a„ is the angle determined in (14.71) . x is the constant cartesian coordinate of the 
unperturbed D6-brane and the only non-zero component of the gauge field / is the one 
written in (14.101) . The equations of motion for the fluctuations can be written in terms of 
the the open string metric (see appendix iBl). which takes the form: 



K 



Gmn dC '^^'^ = 4^ [ ^^^^ ds\^s^ + ] 



(7.2) 



where ds\ij^g^ = —r'^dt^ + dr^jr^ and rfs|,i ^ has been written in (I4.12p and (I4.14p . In terms 
of the metric Qmni the total lagrangian density for the fluctuations takes the form: 



27 F 



sm ar, 



1 3 

_l_ _ nm-p nnq f f _ _ 
~r ^ y y J pq J mn 



uk e/, 



Or 



2 sin a„ sin 



(7.3) 



Let us now write the equations of motion derived from this lagrangian density for the different 
fluctuations. First of all, we write the one corresponding to the fluctuations x of the cartesian 
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coordinates, namely: 



. 



The equations for the scalar ^ and the gauge field are coupled. The equation for ^ is: 



4fc ^9 



3 , 12fc for 



sm a 



i?3 



sm Ur. 



(7.4) 



(7.5) 



whereas the equation for the gauge field is: 



\f9 



d„ 



+ 



I2k 



sin^ ar. 



(7.6) 



7.1 Fluctuation of the cartesian coordinates 



Let us consider in detail the fiuctuations x of the cartesian coordinates. By using the 
expression of in ( IB. 121) . eq. (17. 4p becomes: 



dr{r^ drx) - dlx + r'^ sin^ a„ V|i,i x = 
Let us separate variables in this equation as follows: 

e^^*F,„,,(T^'i)0(r) , 



X 



(7.7) 



(7i 



where Yi^^i^{T^'^) are the scalar harmonic functions for the squashed T^'^ (see appendix [C|) . 
which depend, among other quantum numbers, on two integers li and l2- As shown in 
appendix [Cl the eigenvalues of the laplacian in a squashed T^'^ space with general squashing 
factors depend on /i, I2 and on the R-charge quantum number r (see ( IC.SP ). Remarkably, 
for the particular squashing of (I4.14p the dependence on the R-charge r drops (see (IC.12I) ). 
Actually, by using the eigenvalues of the laplacian derived in the appendix [Cl we get: 



sin^ an VL,i F/ 



{k{h + 2)v + /2(/2 + 2)(l-z/)]r;,, 



(7.9) 



with V being the filling fraction {v = n/N). Using this result, the equation for the radial 
function (h(r) becomes: 



[h{h + 2)u + Z2(/2 + 2)(l-z/)]0 = . 



(7.10) 



Let us now study the behavior of the solutions of this equation for large r. We assume 
that ~ r"^ for large r. By keeping only the leading terms we get that a must satisfy the 
following quadratic equation: 



a{a + 3) = /i(/i + 2)u + hih + 2) (1 - z/) 



(7.11) 
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which has two solutions, namely: 



«± = --±'y/(-) + h{h + 2)u + l2{l2 + 2){l-u) . (7.12) 

Thus, at large r, the fluctuation becomes: 

(f){r) ~ c+r"+ + c_r"- . (7.13) 
In general, in AdS2, if a fluctuation behaves in the UV as: 

X ~ ci r-'"i + C2 r-'"^ , 02 > ai , (7.14) 
then, the dimension of the operator dual to the normalizable mode is: 

1 

A = - + a2-ai . (7.15) 



In our case: 



«i = > «2 = , (7.16) 



and, therefore, the conformal dimension A is just: 



A = ^ + -^(^)' + h{h + 1)v + Z2(Z2 + 2) (1 - I/) . (7.17) 

Notice that A = 2 for l\ = I2 = ^ and, that in general the dimensions are not rational 
numbers. This is, actually, what happens for the supergravity backgrounds that contain the 
T^'^ space [i.e. the conifold-like theories). Moreover, the dimensions depend on the filling 
fraction v. Only in the symmetric representations in which li = I2 = I this dependence on v 
disappears. Indeed, in this case one gets the simple formula: 



A = ^ + + /(/ + 2) , {h^k^ I). (7.18) 

7.2 Coupled modes 

Let us now consider a solution of the equations of motion of the fiuctuations in which the 
gauge field components and Oj are non-zero and given by the ansatz: 

ar = e^^*y,,,,,(fi'i)0(r) , a, = e^^* V,Y;,,,,(f i'^)0(r) . (7.19) 

In addition, we will adopt the ansatz in which the scalar field ^ is also fiuctuating and given 
by: 

e = e^'^*l^„jf^'i);.(r). (7.20) 
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For this ansatz of the gauge potential, the non-vanishing components of the gauge field 
strength f = da are: 

for = lEar, /o. = tEa,, = e*^* {dr4> - 0) V. r/„,,(T^'^) . (7.21) 

We can now write the different components of the coupled equations of motion fl7.5p and 
(17. 6p . The equation of motion (17. 6p for the gauge field for the index n = is: 

3i 

^h,i2^ ^ ^^^^^ ~ ^ Jk^^ sina„a^z . (7.22) 
where c'{_^ i^ is the eigenvalue of sin^a„ V|,i ^ '^hh (see (17. 9p ). namely: 

cr,^,^ = h{h + 2)v + hih + 2) (1 - z/) . (7.23) 
For the index n = r the gauge field equation (I7.6p gives: 

+ Ch,i2 (^r<P - <p) - —r smaniEz = Q . (7.24) 

Finally, for = z we get the following relation between and 0: 

^2 = -r'^dr [r^ (5,0 - 0)] . (7.25) 

One can verify that this last equation (I7.25P is a consequence of (17.220 and (I7.24p . Moreover, 
we can use (I7.22P to eliminate from the equations of motion of the gauge field. The 
resulting equation for is: 

dr{r''dr(t>) + ^0 - <,/,0 - §5 sina, [9, (r^^, z) + ^ z] = . (7.26) 
The equation of motion (17. 5p for the scalar ^ becomes: 

7712 O • 771 At, 

drir'drz) + —z + {3- cl,Jz + -. — = . (7.27) 

Let us now decouple the two equations (17.261) and (17.271) . First of all, we define the differential 
operator O as the one that acts on an arbitrary function ip as: 

772 

0^ = drir^dr^) + ^^Ij . (7.28) 

Next, we define the functions z and t] as: 

1 . ^ 3i . 

z = -— — smanZ , 1] = (p - — — smanZ . (7.29) 
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Then, one can show that fl7.26p and (17.271) can be written as the following matrix equation: 



z 



O- M]\ ^ ) = , (7.30) 



where is the following 2x2 matrix: 

/ cL^,+6 -3 \ 

M = 



(7.31) 



In what follows we will distinguish the case in which ^ from the case in which it 

vanishes and the fluctuation is an s-wave in the internal T^'^ space with h = h = 0. Actually, 
it follows from (I7.22p and (I7.24p that rj = when Ci_^ I2 ~ ^ and, therefore, in this case there 
is only one type of mode, namely that corresponding to z. Notice that (I7.30p is consistent 
with having rj = when Ci_^ I2 — ^■ 

Let us consider from now on the case c^^ I2 ^- this case we have to deal with the 
full matrix equation (I7.30p . The matrix Ai written above is diagonalizable and its two 
eigenvalues are: 

^11 = +3(1 i^/TT^) . (7.32) 



The corresponding eigenfunctions are the following combinations of z and rj: 

^± = i^Jl + cli^ ±l)zTV- (7.33) 
Indeed, one can check that ip± satisfy the following decoupled equations: 

dr {r'dr V^±) + ( ^ - A;^)^ ) ^± = . (7.34) 

In order to obtain the conformal dimensions associated to the normalizable solutions of 
f l7.34p . let us study the behavior of the functions i/j±{r) for large r. By considering functions 
of the type 4'±{r) ~ r", one gets two possible solutions for a and, thus, ip±{r) behaves for 
large r as: 

7/'±(r) ~ + csr-^-^^*' , (7.35) 

where ci and C2 are constants and the exponents a^^^ and 02"^^ are given by: 



X - v' 1 + 4 a1^/ 1 + a/i + 4 Ap/ 
) = ^ , = ^ . (7.36) 

The conformal dimensions A± of the operators dual to the normalizable modes of ip±{r) are 
immediately obtained from the general equation f l7.15p : 

A± = ^ (1 + V'l + 4Ag) . (7.37) 
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Remarkably, 1 + 4 is a perfect square. Indeed, by using fl7.32p one can straightforwardly 
prove that: 

1 + 4A!,1 = (2yi+^±3)'. (7.38) 
Therefore, the conformal dimensions A± can be simply be written as: 

^+ = 2 + , A_ = yr+^ - 1 . (7.39) 

Recall that, to obtain fl7.39p . we have assumed that > 0, which ensures that A_ > 0. 
When c^"^ vanishes then rj = and the matrix equation fl7.30p reduces to the simple equation 
(0 — 6)^ = 0. It is then straightforward to demonstrate that the normalizable mode 
has conformal dimension A = 3, which is the same value as the one that is obtained by 
substituting c^^^^^ = in the equation giving A+ in (I7.39p . 

The dimensions fl7.39p are, in general, not rational and they depend on the filling fraction u. 
However, for symmetric representations in which li = I2 = I the dependence on v disappears 
and one has: 

cl = 1(1 + 2) =^ ^/l + ^i = l + l. (7.40) 
Therefore, in this symmetric case the dimensions are rational and are given by: 

A+ = 3+/, (/>0), A_ = / , (/>0), (7.41) 

where we have already taken into account the / = case. 



7.3 Internal gauge field modes 

We now consider fluctuation modes of the gauge field in which the only components which 
are non-vanishing are those along the directions of the internal T^'^ space (they will be 
denoted by at). In this case the non-zero components of the field strength / are: 

foi = dotti , fri = drtti , fij = di aj - dj ai . (7.42) 

The equations of motion (17. 60 for the gauge field for n = 0, r are just: 

9o(V'ai) = , dr{W'ai) = , (7.43) 

where: 

V^a, = (Vl^^^ a,) . (7.44) 
Thus, clearly we should require the fulfillment of the following transversality condition: 

V^Oi = . (7.45) 
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We now write the equation (17. 6p for n = i, namely: 



dr [ dr a, ] - dl a, 1 + ^J^^ (^fc^i* _ g'a'' ) ] = . (7.46) 



The last term in this equation can be rewritten as: 

^ dk [v^ {d^a' - d'a^ ) 1 = V^V^ a' - R\ a!" - V Vfc , (7.47) 

where is the Ricci tensor of the T^'^ space. Taking (I7.45P into account, the last term in 
(I7.47P vanishes and the right-hand side of (I7.47P can be written in terms of the Hodge-de 
Rham operator Ai, which acts on a vector field with components fi as: 

Ai/, = -VfcVV. + R-fk . (7.48) 
Thus, the equation of motion (17.460 becomes: 

dr[r^drai] - ^ dl at - 4z/ (1 - z/) Ai = , (7.49) 

where we have taken into account that sin^an = 4z/(l — z/). Let us next separate variables 
in tti as: 

a, = e^^*y,(Ti'i)0(r) , (7.50) 

where Yi{T^'^) is a vector spherical harmonic for the T^'^ space which satisfies the transver- 
sality condition: 

V*ri(fi'i) = . (7.51) 

Let us diagonalize the Hodge-de Rham operator Ai in the space of the T^'^ vector harmonics 
and write: 

4z/(l-z/)Aiy, = Ar, , (7.52) 

where the eigenvalue A will depend on the quantum numbers li and Z2, as well as on the 
filling fraction 1/ (see below). By plugging the ansatz fl7.50p on the fluctuation equation 
( I7.49p . one gets the following equation for the radial function 0(r): 

(94r2a,(/)] + — - A(/) = . (7.53) 

This equation can be solved for large r by taking ~ r", with a being some exponent. 
Taking into account that the term with the energy is subleading at large r, one gets that the 
exponent a must satisfy the following quadratic equation: 

a (a + 1 ) = A , (7.54) 
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which has two solutions, namely: 



-l±v/rT4A 

a± = ^ • (7.55) 

Thus, 0(r) behaves as f l7.14p with ai = — a+/2 and 02 = — a_/2, z.e. as in f l7.16p . Therefore, 
the conformal dimension of the operator dual to these fluctuations is: 

A = ^ (^1 + v^Tma) . (7.56) 

Let us analyze the different dimensions obtained for the various vector harmonics for the 
T^'^ space that were studied in appendix [Dl As shown in this appendix, the results depend 
on whether the filling fraction u is arbitrary or z/ = 1/2. These two cases will be studied 
separately in what follows. 

7.3.1 Arbitrary filling 

As argued in appendix [D] in this case we should take vanishing R-charge r and the eigenvalues 
A depend on the remaining two quantum numbers /i and 12- Actually, there are two series 
of eigenvalues, which were denoted by Ai and A2 in the appendix. Let us consider the first 
of these modes and let us denote by Ai the corresponding value of A in (I7.52p . From the 
value of Ai written in (ID. 281) . we get: 

Ai = . (7.57) 

The associated conformal dimension Ai is given by: 

1 



which is, in general, not rational. Moreover, only for li = I2 = I the dimension Ai is 
independent of the filling fraction since dd = l{l + 1). 

The second class of eigenvalues in this arbitrary filling case, denoted by A2'^^ in (ID.28p . 
lead to a value of A given by: 

a(±) = 1 + ± + 4z/3 + 2) + 4(1 - z/)3 /2(/2 + 2) . (7.59) 

The corresponding conformal dimension is: 



Ai^^ = 1 + ^1 + cl,^ ± V^ + AuH,{h + 2) + Ail-i.rhik + 2) . (7.60) 
We will see below that this formula simplifies greatly if u = 1/2. 
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7.3.2 Half filling 

When V = 1/2, it was shown in appendix |D] that there are transverse modes for r 7^ 0. As 
in the generic u case, there are two types of modes. The eigenvalues of the first type were 
denoted by A 1 in (]D.29P and they lead to the following eigenvalue in (17.521) : 

= ci,,i, ± r , (7.61) 

with Qj^/j being the value of cj"^ ^ = V^; ^-e.: 

- _ u=i/2 h{h + 2) , h{h + 2) 

ci^M = = ^ + ^ • (7.62) 

The dimension of these modes are non-rational and given by: 

^^^^ = 1 + y\ + ^h,h±r. (7.63) 



For r = these dimensions reduce to the ones written in fl7.58p for z/ = 1/2. 

Next, let us consider the modes with 
fl7.52p is independent of r and given by: 



Next, let us consider the modes with eigenvalues for which the eigenvalue defined in 



Af ) = 1 + Qi,i2 ± + ^h,h ■ (7-64) 

Remarkably, one can verify that 1 + 4:A.^2^^ can be written as a square. Indeed, one can 
straightforwardly check that: 

1 + 4A^^) = (2v/l + c\,,,±l)' . (7.65) 

Thus, the associated conformal dimensions Ag^'' and A2 ^ are just: 

A^+) = 1 + y/1 + Q,,,, , A^-) = + • (7.66) 

These dimensions are generically not rational. However, when h = I2 = I, one has that 
a/c; I = + 2) and, thus, 1^1 + q^i = / + 1. Therefore, in this h = I2 = I case the dimensions 
associated to these modes are rational and given by: 

a5+^ = 2 + / , A^"^ = 1 + / , (/i =l2 = l) . (7.67) 

8 Dimer fluctuations 

Let us now consider the fluctuations around a connected conflguration of the type studied 
in section El We will restrict ourselves to analyze the case in which the background is 
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supersymmetric, i.e. for vanishing temperature. Recall that in this case the D6-branes do 
not reach the origin at r = and they extend along the Minkowski directions (say, along 
the coordinate x). For this reason we will have to distinguish between longitudinal and 
transverse fluctuations in the {x, y) plane. By computing the lagrangian for generic quadratic 
fluctuations one realizes that most of the modes are coupled and the diagonalization of the 
equations of motion is rather involved. However, one can verify that the fluctuation of 
the transverse position in the Minkowski space is decoupled from the other modes. In this 
paper we will limit ourselves to studying these transverse fluctuations. The corresponding 
lagrangian density and equation of motion are derived in appendix |B] (eqs. f lB.37|) and 
( IB. 381) . respectively). More explicitly, the equation of motion for the transverse modes can 
be written as: 

dr{r' Vr'-rtdrx) " r^^ d^oX + r^-^ sin^ a, V|,,i x = , (8.1) 

where is the minimal value of the r coordinate attained by the unperturbed configuration 
of the brane. Notice also that any value of r in the interval [r^,, +oo) is reached twice by the 
brane. 

Let us separate variables in (18. ip as in (17. 8p and let us define the rescaled radial coordinate 
p and the rescaled energy E as: 

r - E 

p = - , E = - . (8.2) 

In terms of these quantities, the equation of motion for the fluctuation becomes: 

P' u P' 



9p(/y7^5p0) + ^'^^0 - = , (8.3) 

VP - 1 VP - 1 

where c'{_^ i^ has been defined in (I7.23p . 

When Ci_^ I2 ~ ^ (*^^ equivalently for the s-wave modes with li = I2 = 0), the fluctuation 
equation (18. 3p can be exactly solved by means of a change of variables found in |2ll 125] . 
Indeed, let us introduce a new variable 77, related to p by means of the equation: 

V = r . (8.4) 

By this change of variables p G [1, +00] is mapped into r] G [0,f]{E) ], where fi{E) is defined 
as: 

viE) = r -"^^ , . (8.5) 

Next, we introduce a new function as follows: 

$(r^) = P^ <p{p). (8.6) 
VP^ + E^ 
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By explicit calculation one can prove that the equation satisfied by $(77) is simply: 

(P^ 1 



+ -E'{E^-1)^ = , (8.7) 



which can be immediately integrated in terms of trigonometric functions. By imposing that 
$ vanishes at 77 = 0, f]{E), we get the following quantization condition on the energy E: 



E - 1 f]{E) = rnr , n = l,2, (8.8) 

This equation must be solved numerically in order to get the energy levels. 

Let us come back to the analysis of the equation fl8.3p for general values of the Kaluza- 
Klein quantum numbers h and l2- We will show that this equation can be mapped to the 
Schrodinger equation of a quantum integrable model, namely the so-called Inozemtsev BCi 
model. The first step to prove this result is changing the radial variable p by a new variable 
z, defined as: 

z = . (8.9) 

Then, one can check that the fiuctuation equation (18.31) can be written as: 



d'^ _^ 1 



dz^ 2 



1 3 1 



z z — 1 z — 2 



2 



1 cY , (z-l) + E 

d^(t)--Jlhl > ^0 = 0. (8.10) 

^ 4 z{z-l){z-2) ^ ^ ' 



This equation is a Heun's equation, whose more general expression is given in appendix 
lEl (eq. flE.ip ). As reviewed in appendix |El the Heun equation is an ordinary differen- 
tial equation with four regular singular points whose location and characteristic exponents 
are parametrized by several numbers, which we denoted by a, /3, 7, 5, e, q and d (see 
(]E.1|) ). The solution of the general Heun equation ( lE.ip which is regular at z = is denoted 
El{d,q\a,p>,-i,b\ z). 

By direct comparison of ( 18.10P and (IE. II) . one gets that in our case d = 2 and that the 
coefficients 7, 5, e and q are: 

7 = -, = - , e = - q = ■ . 8.11 

' 2 ' 2 ' 2 ^44 ^ ' 

Notice also that the remaining parameters a and /3 should satisfy (see flE.2p ): 

a + P = + 6 + e - 1 = I , a/3=-^. (8.12) 

In order to determine a and /3 we use the fact that that they are the roots of the quadratic 
equation: 

3 c'^ 

- {a + 13) X + al3 = - -X - = , (8.13) 
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which are: 

a-^. e-'-^. (8.14) 

with A being the conformal dimension (I7.17p . Then, the solution of the fluctuation equation 
which is regular at p = 1 is: 

+ (A + l)(2 - A) A + 12-A13 
0(P) = Hl{2, "--^ (8^15) 

When ^ = we can rewrite this result in terms of a hypergeometric function (see flE.SP ). 
namely: 

^(P)U.„ = ^'(^.^.^i (8.16) 

As reviewed in appendix [El with an additional change of variable, the Heun equation 
can be converted into a Schrodinger equation for a particle moving under the influence of 
a potential which is a combination of elliptic functions. For a generic Heun equation this 
potential was written in eq. flE.lSp . For our particular case of f lS.lOp the hamiltonian takes 
the form: 

H = --^ + A(A-l) p{x) , (8.17) 

where p{x) is the Weierstrass function for periods (2a;i,2a;3) = (see eq. flE.6P for its 

general deflnition). The relation between the variable x in fl8.17p and the coordinate z of 
the Heun equation has been written in ( ]E.25p . It is interesting to write here the relation 
between our original radial variable r and the new variable x. One gets: 

Lr 



where L is the length of the dimer (see eq. ( IB.32P ). By means of this change of variable the 
equation for the fluctuation can be written as: 

H<p = 8<p , (8.19) 

where is the same function as in f l8.10p but now considered as a function of x. Moreover, 
the energy £ of the Schrodinger problem f l8.19p is related to the energy E. Indeed, by using 
dESI]) one gets: 

= ^ 8 . (8.20) 

.r(i)_ 

Notice that the Schrodinger equation f l8.19p for the hamiltonian f l8.17p is nothing but the 
Lame equation. When A G Z the eigenfunctions and band structure of the Lame equation 
have been obtained long ago [26]. Notice, however, that in our case A G Z only when 
/i = ^2 = (see (EnD). 



29 



The coordinate x in f l8.17p takes values in the interval [0, 1], where the Weierstrass function 
'p{x) with periods (2a;i,2a;3) = (1,^) is real and positive and has a minimum at x = 1/2. 
Moreover, p{x) — )■ +oo when x — 0, 1, which means that the eigenvalue problem has a 
discrete spectrum consisting of an infinite tower of eigenvalues. Although many results are 
known for these type of integrable models, we will leave this analysis for a future work and 
here we will content ourselves with estimating the eigenvalues by means of the WKB method. 
Actually, this estimate can be performed directly in the variables of eq. flS.ip . Indeed, by 
using the equations written in p7], we get that the energy levels can be estimated as: 



c 



(n + 1) (n+ A) , 



where C, is given by the following integral: 

dp 



c 



1 



4V27r 



{n > 0) 



(8.21: 



1 1 2 



(8.221 



Let us write this result in terms of the length L of the dimer. The relation between L and 
was written in eq. ( ]B.32[) . By using this result, one finds that the un-rescaled energies En 
are given by: 



E ~ 



[r(i)]'i 



(n + l) (n + A' 



(n > 0) . 



(8.23) 



9 Summary and conclusions 

Let us briefiy recapitulate the main content of this work. We have studied the addition of 
localized fermionic impurities to the ABJM Chern-Simons-matter model in 2+1 dimensions. 
In the holographic approach the impurities are added by means of D6-branes extended 
along the holographic coordinate and wrapping a squashed T^'^ space at a fixed point of 
the Minkowski spacetime. The coupling of the D6-brane to the RR field of the background 
induces an electric worldvolume gauge field which prevents the collapse of the wrapped brane 
and must be quantized accordingly. 

The number n of defect fermions fixes the size and the internal deformation of the T^'^ 
space. We have studied the thermodynamics of the system by analyzing the D6-brane 
action in a black hole background and we have determined the entropy and free energy 
of the impurity. We have also studied connected configurations in which the branes are 
hanging from the boundary and we have shown that these configurations dominate over the 
unconnected ones at low temperature, giving rise to a dimerized phase. 

A crucial point in all the approach is the stability of the D6-brane embeddings used 
to introduce the holographic impurities. In order to study this aspect of the system, we 
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have performed a detailed analysis of the fluctuations of the D6-brane probe in the ABJM 
background and we have determined the spectrum of conformal dimensions of the dual 
operators in the defect theory. This required a detailed analysis of the KK modes of the 
probe brane and of the harmonic analysis on the squashed T^'^ space which uncovered a very 
rich structure. We also studied a subset of fluctuations of the connected configurations and 
pointed out a connection with an elliptic quantum integrable model. 

There are many unanswered questions raised by our results which could be the subject 
of future research. First of all, to complete the holographic setup we should have a more 
detailed description of the field theory side of the correspondence. In particular, we would 
like to determine the concrete action of the fermionic fields of the impurities and their 
coupling to the different bulk fields. This would allow us to write the precise dictionary 
between fluctuations of the probe and impurity operators. Moreover, in order to complete 
the spectroscopic analysis of section [7] we would have to study the fermionic fluctuations 
of the D6-brane probe and we should verify how the different modes are accommodated in 
supermultiplets. 

All our results have been obtained in the probe approximation. As argued in [10], there 
are some physical effects which are not captured by this approximation and that require 
taking into account the backreaction from the impurity. In this sense it would be interesting 
to find a generalization of the bubbling geometries [28] for the ABJM Wilson lines, similar 
to the ones found in [291 EQl [31] for the AdS^ x geometry. 

On general grounds it is also of great interest the extension of our results to less super- 
symmetric configurations and/or to systems in which the ambient theory is not conformal. 
In the first case one could study impurity theories constructed from the brane embeddings 
of [32]. In the case of non-conformal bulk theories we have studied the addition of impurities 
to backgrounds generated by Dp-branes, with p ^ 3. The results will be presented elsewhere 
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A BPS configurations 



The Euler-Lagrange equation of motion derived from the hamiltonian density f l3.20p is: 



d_ 

dr 



\j ( sina)*" + C„i 



{a'Y 

^=^= (sina)^ AnCal 



sin a) + Cn(a)^ 

(A.l) 

where A„(q;) is the function defined in (14.61) . In the main text we have studied solutions of 
these equations with constant a. We will now find a more general class of solutions of f lA.ip . 
Instead of dealing with the Euler-Lagrange equations we will follow a different strategy by 
establishing a BPS bound for the energy (see [371[38]). The configurations that saturate this 
bound solve the equations of motion and will be characterized by a first-order differential 
equation which can be integrated in analytic form. We will start by rewriting the hamiltonian 
density fl3:20|) as: 

where Z and y are given by: 

Z = r a' [(sina)^ cos a — C„(a;) sin a] + (sina)^ + C„(q;) cos a , 



y = (sina)^ cos a — C„(a) sin a — r a' [(sina)^ + C„(a) cos a] , (A. 3) 

and Cn{<y) is the function defined in (I3.14p . Clearly, the energy functional H satisfies the 
bound: 



which is saturated by the configurations that satisfy 3^ = 0. Interestingly, for any function 
a(r), the function Z can be written as a total derivative: 



dr 



sm a 



Ar,(a) cos a 



(A.5) 



As a consequence of (1A.5I) . only the boundary values of a{r) matter when one evaluates the 
right-hand side of f lA.4p and, thus, the functions a{r) which saturate the bound correspond 
to D6-brane embeddings which, for given boundary conditions, minimize the energy. Let 
us now analyze in detail these minimal energy configurations. In order to write the BPS 
condition 3^ = in a suitable way, let us recast y in the form: 

y = sinaA„(«) — r a' [(sina)^ — A„(a) cos a] . (A. 6) 

Using this result, the BPS condition can be written as: 

, 1 sinaA„(a) 



a 



r (sina)^ — A„(a) cos a 
32 



(A.7) 



Several observations concerning (1A.7P are in order. First of all, notice that the solutions 
of (I A. 71) with constant a are those with a = 0, tt or with a being one of the zeroes of the 
function A„(a;). Thus, our flux tube embeddings are certainly a particular solution of the 
BPS equation flA.7p . Moreover, one can verify by direct calculation that any solution a(r) 
of (1A.7I) also solves the second-order differential equation of motion ( lA.ll) . To prove this fact 
it is quite useful to demonstrate first that: 



, (sina) + Cn(a)2 

1 + [a'Y = , . (A. 

' BPS (smaj^ — An cos a 



Furthermore, by using flA.8|) one can also prove that the electric field for a BPS configuration 
is given by: 

\bps = ^ '^°^") • (^•^) 

Actually, one can easily demonstrate by using f l3.16p that the first-order BPS equation ( 1A.7I) 
for the embedding is equivalent to having an electric field given by (IA.9p . Moreover, although 



we have not verified it, it is likely that the first-order equation flA.7p and the electric field 



( 1A.9P could by derived from the kappa symmetry of the probe brane, in a calculation similar 
to the one done in [39] for the A/" = 4 baryon vertex. 

By using the fact that dAn{a)/da = 2 sin a, one can easily show that the BPS equation 
can be integrated exactly in the form: 



r = Cy^^^, (A.IO) 
sin a 

where C is a constant of integration. By redefining C and using the expression of A„(a;) in 
(14. 6p . this solution can be written as: 



cos a + — 1 



2n 

C ^ ^ . (A.ll) 

sin a 

Notice that, in the solution flA.lip . the coordinate r — )■ oo when a — )■ 0, which corresponds 



to spike of the D6-brane probe which reaches the boundary of AdS^. Moreover, the angle a 
takes values in the interval < a < a„, where the a„ are the angles defined in (14. 5 p (when 
< n < N the brane reaches the origin r = when a = a„) . By looking carefully at the 
energy of the spike one can check that it corresponds to n fundamental strings reaching the 
boundary of the Anti-de-Sitter space. 

In the particular case n = N the maximum value of a is a = tt (z.e. < a < vr in this 
case) and the D6-brane wraps completely the CP^. This configuration corresponds to the 
baryon vertex of the ABJM model. In this n = N case the function r{a) can be simply 
written as: 

K«) = , < a < vr , (A.12) 

sm (2) 
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where Tq = r{a = vr). 



B Fluctuations of the probe D6-branes 

In this appendix we analyze the small perturbations around the flux tube and dimer con- 
figurations at zero-temperature. The goal is to find the second order lagrangian for these 
fiuctuations that was the starting point of sections [7] and [HI We start by studying the 
fiuctuations of the impurity configuration of section HI 



B.l Impurity fluctuations 



Let us perturb the configuration of the D6-brane probe as in (17. ip and let us expand the 
DBI+WZ action to second order in the perturbations ^, / and x- We begin by obtaining 
the different components of the induced metric gmn at second order in the fiuctuations. We 
write: 



9 = 9 + 9 



[B.l] 



where g is the zero-order induced metric written in (14. lip . When any of the indices m and 
n is outside the T^'^, the metric perturbation Qmn takes the form: 



9mn - 



dmi dni + dmX dnX 



(B.2) 



On the other hand, if both indices belong to the T ' , one has to expand the sin | and cos | 
factors appearing in the internal metric up to second order in ^ and one gets: 



9ij 



Ak 



d^id.i + r'dad.x + ig)'' +^91 



(B.3) 



where g'^^'> and g'^'^^ are given by: 



gff d^^ 



gff d^^ 



— sin ( dOf + sin^ 9i dip\ ) + sin a„ ( dO\ + sin^ 6*2 dip\ ) -|- 

-|- sin an cos an{dx + cos 6i dipi + cos 62 d(p2 ) ' 

— cos an ( dOl + sin^ 6*1 dipl ) + cos a„ ( dO^ + sin^ 62 d(p\ ) + 

+ (sin^a;„ — cos^ a„) i^dx + cos 9 1 dip i + cos 6^2 (iv92 )^ . (B.4) 
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In order to expand the DBI D6-brane action (13. 3p . we notice that the Born-Infeld determinant 
can be written as: 



V- det(^ + J^) = ^- det ( ^ + / ) Vdet(l + X) , (B.5) 
where the matrix X is given by: 

X ^ {-g + f)-' {g + f) . (B.6) 
To evaluate the right-hand side of eq. (IB. 51) . we shall use the expansion: 

Vdet(l + X) = 1 + ^TrX - ]ttX^ + I (TrX)' + o{X^) . (B.7) 

2 4 o 

Moreover, in the inverse matrix [g + f) ^ we will separate the symmetric and antisymmetric 
parts: 

[g + /)"' = G-' + J , (B.8) 

where J7 is the antisymmetric component. The symmetric matrix Q is the so-called open 
string metric and is the one that naturally shows up in the fluctuations of the worldvolume 
when gauge fields are turned on. 

The matrix g + f has a block structure. By computing the inverse in the Or sector {i.e. 
in the AdS2 part of the geometry), one gets: 

It follows from (IB.9P that the non-vanishing elements of J are: 

j- = -:r- = ^^. (B.io) 

it"^ sma„ 

In terms of the filling fraction v = n/N, the above expression of J7'°'" can be written as: 

0. ^ _jro ^ k N{N-2n) ^ k 1 - 2u 

W n{N -n) ^ ' ' 

Moreover, the non-vanishing elements of the inverse open string metric are: 

r-OO ^ rrr ^k ^k ~ij m i 

and, thus, the open string metric takes the form: 
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where ds\^g^ = —r'^dt^ + dr^jr^ and ^ has been written in (14.12^ and (I4.14p . 

Let us now compute the different terms on the right-hand side of flB.7p . The matrix 
elements that contribute to Tr(X) are: 

4/c cotan „ 1 1 (r. \i , cota„ 

R-^ sm an sm a„ sm a„ sm a„ 

4/ccota„ 2 , . .2 cota„ 

tv" sm an sm a„ sm a„ sm a„ 

4A- 

= + + (-M^'^);. + ^^'^..e^.e + , (B.14) 

where j, are indices along the internal directions and the matrices ^A^°^ are defined as: 

~.fc^(a) ^ (a =1,2), (B.15) 

with g being the T^'i metric fH:T2|) and (f""^ are the metrics written in flB.4p . One can verify 
that the matrices MS""^ are, in fact, diagonal. Actually, ordering the directions of the internal 
manifold as: 

y = (^i,^2,¥^i,¥^2,x) , (B.16) 
it is easy to check that A^^^^ and A^^^^ are just: 

MS^^ = diag - tan (^^^ , cot (^^^ , - tan (^^^ , cot (^^^ , 2 cot a, 

(2) ,. ( COSa„ COSa„ COSa„ COSa„ o i\ it:, ^n\ 

-^"H- 4^^' 4^1^' -4^^' 4^^ ""-^j- (^-^'^ 
It follows that the traces of the A^'s are: 

Tr (AlW) = 6cota„ , Tr (Al^^)^ ^ 3pQ^2 _ 1 . (b.18) 

From these expressions it is now straightforward to compute the trace of X, with the result: 

Tr(X) = -.^^^^^/o, + 6cota„e +(3cot2a„ - 1) + 
sma„ 

+ + ^r2e;-"9„x5nX- (B.19) 

4A; 4k 

In order to calculate the trace of X^, one needs to compute the non-diagonal elements of X. 
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To the order we are working, it is enough to calculate them at first order. One gets: 
4fc 1 „ Ak 

sin^ ar 



" i?3 r^sin^a^ ' " --^^ " ^' 

_ 4fe 1 ^ 4A;cota^ _ 



Moreover, by using that: 

Tr (M« A^(i)) = 4 (3 cot^ a„ + 1 ) , (B.21) 

one gets: 

Tr(^^) = ^ '""''""^ V o^. + 4(3cot^a„ + 1 ) + 
sm a„ 

I / 2 2r^ X 

/t° V r"^ sm sm a„ / 

With these results one can now readily compute A/det(l + X), namely: 

4A; cotQ;„ . ^ ^ i?^ 1 



Vdet(l + X) = 1 - —-. ^/o. + Scota^e + 77 77 e^'"" ^ 5„ ^ + 

it'^ sma„ 4A: 2 

+ ^(2cot2a„ -l)e - ^^^e/o. . (B.23) 
2 it'^ sm On 

Let us now use this result to compute the DBI action for the fluctuations. First, we notice 
that: 



e-^^J-det{g + f) = ^ sin«„ , (B.24) 

where g is the determinant of the metric of the squashed T^'^ which, according to (I4.12p . is 
given by: 

= -— sin^i sin ^2 • (B.25) 
8 

Therefore, if we represent the DBI action as: 

Sdbi = j Cdbi d\ , (B.26) 
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the DBI lagrangian density is given by: 

C-DBi = -Tq sina„A/^ ^/det{l + X) , (B.27) 

where A/det(l + X) has been written in terms of the fluctuations in (]B.23p . 

The WZ lagrangian density has been written in (13.81) . Let us write in this equation 



cos On + for and let us expand the function C{a) as: 



9 

C{a) = C(a„) — 3(sina„)^^ — - (sin a^)^ cos a„ . (B.28) 



Then, one gets: 



^wz = -Tq sin 6^1 sin 



^C(a„) cosa„ + C(a„) /or 



37^3 g^3 -, 

— — sin^a„ cosa„^ — sin^ a„ cos^a„^^ - 3sin^a„^/or • (B.29) 

The total lagrangian density for the fluctuations is just the sum of ( IB. 271) and ( IB. 291) . After 
neglecting terms that are constant or total derivatives, one can easily prove that the total 
lagrangian density is given by (17. 3p . 



B.2 Dimer fluctuations 

Let us now study the fluctuations around the connected configuration given by the ansatz 
(16. ip in the zero-temperature ABJM background of section [2l First of all, we notice that the 
function x{r) is given by: 

? = ± 4=, (B.30) 

with being the minimum of the coordinate r reached by the brane. Eq. flB.30P is immedi- 
ately obtained from (I6.10p by taking /(r) = 1 and by redefining the constant A as = 
(see (I6.14P ). The length L of the dimer is now: 

L = 2rl -^j== . (B.31) 

By performing explicitly the integral on the right-hand-side of ( ]B.3ip . we get: 

2y27rf 1 



(B.32) 
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The induced metric in this case becomes: 
where h{r) is the function: 



4k 



dt"^ + dr'^ + dsL i 



h{r) 



1 - 



(B.33) 



(B.34) 



Notice that in this case the metric is of the form AdS2 x T^'^ only asymptotically in the 
UV. The worldvolume gauge field strength for this hanging configuration depends on the 
holographic variable r. Indeed, from fl6.13p one easily gets: 



for 



Ak 



r cosa^ 



(B.35) 



(Compare this result with (14.101) ). Proceeding as in the case of the unconnected configura- 
tion, it is straightforward to obtain the open string metric Qmn for this case. One gets: 



m n 
QmndCdC = 



sin^ ar. 



^-r'dt' + ^ dr' ) + ^4,, 



(B.36) 



The lagrangian density for a generic fluctuation of the dimer is highly coupled and, for this 
reason, we will not try to study here the dimer fluctuations in full generality. In this paper 
we will restrict ourselves to a set of fluctuation modes which are decoupled from the others. 
In these modes the only coordinate that fluctuates is the cartesian coordinate transverse 
to the dimer direction in the Minkowski space {i.e. y = Ho + Xi with hq constant). It is 
straightforward to expand the DBI+WZ action of the probe and to obtain the lagrangian 
density to second order in x- Oi^s gets: 



i?9 



27 P 



sm Ur 



Ak 2 



9m X X 



(B.37) 



By comparing this expression with the first term in (17.31) we notice that in ( IB.37P there is an 
extra factor of y/h, and the difference between the lagrangians is similar to the one between 
the open string metrics (]B.36P and (]B.13p . The equation of motion derived from (1B.37P is: 



X 



(B.38) 



This equation is studied in detail in section [81 
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C Laplacian for the squashed T^'^ 

Let us consider a five- dimensional metric of tlie type: 



A{d9l + sin^ 9id(f)l) + B{ ddl + sin^ 9 



C {dx + cos 6i d(j)i + cos 62 



(C.l) 

where A, B and C are constants. Tlie laplacian operator for the metric ( IC.ip acts on a scalar 
function H as: 

(C.2) 



where, for convenience, we have defined the operator Aq = — (Aq is the Hodge-de-Rham 
operator acting on zero- forms). Let us now define the following operators: 



1 



sin 9, 



deXsmOide^) + 



sin 9, 



cot 9i 



1,2) 



Then, the laplacian Aq can be expressed in terms of 2 ^"^^ ^"r 

The laplacian operator for a round three-sphere is given by: 

1 



§3 



sin 6' 



de{sm9d0) + 



1 \ 2 

dd, — cot9dy) + dl 



sin 6* 



(C.3) 



(C.4) 



(C.5) 



and its eigenvalues are — /(/ + 2) with / G Z. It follows that the operators V? are related to 
the laplacian of a three-sphere parametrized by (6*,, 0j, x) by means of the relation: 



52 . 

X 



(C.6) 



The eigenfunctions of d'^ are of the form e*2^ with r G Z and the eigenvalues are of the form 
— Then, the eigenvalues of are: 



kik + 2) r2 
4 4 



Therefore, the eigenvalues of Aq are: 



HQ{li,l2,r) = — h 



4A 



AB 



A B 



1 1 r 



C 



(C.7) 



(C.8) 
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Let us apply this result to the case of the squashed T^'^ metric written in fl4.12p . In this 
case the values of A, B and C are: 



A = cos 



2 



B = sin 



2 



COS 

2 2 



C = sm 



2 ' 2 ' 

In terms of n and the actual values of the coefficients A, B and C are (see ( 14.14p ): 



A 



These values satisfy: 



N -n 
N 



B 



n 
N 



C 



n{N — n) 



1 1 _ 1 
A^B~C' 

and, therefore, the eigenvalues of the laplacian are independent of r and given by: 



Ho- J 



k{h + 2) , 12(12 + 2) 



N -n 



+ 



n 



(C.9) 

(C.IO) 
(C.ll) 

(C.12) 



D Harmonic calculus 

In this appendix we develop the harmonic calculus for a squashed T^'^ space. We will follow 
closely the algebraic approach of refs. |10l HI], in which the T^'^ space is realized as a coset 
of the type SU{2) x SU{2) lUniX)- We denote by Ti, T2 and T3 the generators of the first 
SU{2) and by Ti, T2 and T3 the generators of the second SU{2). The UniX) is generated by: 

= T3 + Tg . (D.l) 

Moreover, let us introduce T5 as: 

T5 = T3 - fa . (D.2) 

In the following we will adopt the notations used in refs. jlQl HI]. The indices a,b,c - ■ ■ = 
1, ■ ■ ■ ,5 will be used to denote a coset direction, i,j,k--- = 1,2 will refer to the directions 
along Ti 2, while the indices r, s, t ■ ■ ■ = 3, 4 will denote the directions along Ti 2. The internal 
indices have a negative definite metric. We will also work with the combinations T± and T±, 
defined as: 

T± = Ti ±iT2 , f ± = f 1 ±1X2 . (D.3) 

In order to describe the coset geometry, we will use a basis of rescaled vielbeins V"" = 
(y^ , V , V^) with the rescaling factors r(a) being given by: 

r{i) = a , r(s) = b , r(5) = c . (D.4) 
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The vielbeins satisfy the torsion-free condition: 

dV = B"^ Vb , (D.5) 

where B"''^ is the Riemann connection one-form [B"''' = —B^'^ for a,b = i,s,5). The 5*0(5) 
covariant derivative is defined as: 

V = d+B'''Tab, (D.6) 

with Tab being the 5*0 (5) generators. As explained in detail in |10l E], the covariant deriva- 
tive V acts algebraically on the coset representatives, which can be taken as the basic har- 
monics. Indeed, let M"^^ denote the part of B""^ orthogonal to the so-called iJ-connection. 
Then, the covariant derivative acting on the coset space can be written as: 

V = -r{a)V''Ta + M^^T.^ . (D.7) 

In order to write explicitly the form of "D, let us notice that any 5*0(5) representation can 
be branched with respect to its U{1)h subgroup. In particular, the vielbein V"", which 
transforms in the 5*0(5) vector representation, can be decomposed into five one- dimensional 
fragments = (y^,V'^), = (V^jV^) and with U{1)h charges given respectively 
by (1,-1), (1,-1) and 0, Then, the components of the covariant derivative V acting on 
harmonics are: 

■ 

Vi = -aTi - — e/ Tgj , 
V, = -hT, + Yc^sTu , 

V, = -cT, - 2(c-^)Ti2 + 2(c-^)t34 . (D.8) 

The second-order Laplace-Beltrami operator is just VaV"-. Let us first compute its eigen- 
values acting on scalar harmonics. In this case the 5*0(5) representation is trivial and T>aT>'^ 
reduces to Aq, where Aq is given by: 

Ao = -a^TiTi - h'^TsTs - c^T^T^ . (D.9) 

In order to find the eigenvalues of Aq, let us consider the action of the T° generators on the 
spherical harmonics. The latter are denoted by Yq^'^'^\ where j and I are the spin quantum 
numbers of the two SU{2) {j, I G Z/2), r is the charge under the t/(l)_R generated by T5 and 
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q is the charge under U{1)h- The generators T±, T± and T5 act on the F's as follows: 

± 



± 



From this equation one gets: 

J", y 



-^s -^s -fg 



(g + r)2 
4 

(g - r)^ 
4 



9 



Then, one can demonstrate easily that: 



(D.IO) 



2 y(i,Z,r) 
5 



2 2 
2 „2 t2\ ''^^ /„2 , 7,2\ 9 / 2 7,2\ 



(D.ii: 
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a'j{j + l) + b'lil + l) + i4c'-a'~b')--{a' + l^')---{a'-f')--^ 

(D.12) 

The scalar harmonics correspond to taking g = 0. In this case the operator Aq acts as: 



{j,l,r) 



HnYr 



(J,l,r) 



-f ' 



where the eigenvalue Hq is given by: 



4p2 _ „2 _vi 

Ho ^ a' jU + l)+b' 1(1 + 1) + J . 



(D.13) 



(D.14) 



Let us now consider the Hodge-de-Rham operator Ai which, in the notation of |40| 141] . 

acts on the vector harmonics Y(a) as: 

Air(,) = V.V'Y^a) + 2i?,'F(b) , (D.15) 

where is the Ricci tensor. Actually, for the SU{2) x SU{2)/U{1) coset metrics, the Ricci 
tensor is block diagonal and its components are given by: 



St 



(D.16) 



2 16c2 ; ' "5 8c2 ■ 

Generically, Ai acts non-diagonally on the vector harmonics Yj^a). Let us represent its action 
in terms of a matrix Xi: 

Air(,) = . (D.17) 
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We will represent the generators of 5*0 (5) in the fundamental representation by means of the 
matrices {Taby^ = | i^a^b ~ ^a^b)- Moreover, in order to write the components of Xi it is 
convenient to use as in [lOl HI] a complex basis for the vector representation, by defining the 
components (±) = 1 ±i2 and (±) = 3 ± z 4. In this basis the only non-vanishing elements 
of the matrix Ai are: 



Ao + + — T2tc[c - — T5 



4c ^ ' 



M 



± 

± 



Ao + + - ±2ic[c - -T5 



2c ^ ' 



b^ ^ 
2c ' 



Ao + 



(g^ + by 
8? 



(D 



The elements of Ai contain generators which act non-trivially on the harmonics Y(a). 
us represent this action in terms of the matrix Ai as: 



Let 



(D.19) 



To compute the values of the different elements of Ai we have to specify the particular 
harmonics which compose our vector Y(a). Let us follow again |10l HI] and define: 



^±1 ' 



Y. 



^±1 



(±) = ^±1 ' ^{±) 

We organize the harmonics as the following vector: 



Y< 



(0) 



Yr 



{j,l,r) 



Y( 



(a) 



Then, then the non-zero elements of Ai^ 

,2 



Yr 



(-) 



Yr 



a 



HoT-r 
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\ ym J 

are: 



c? ( r 



(D.20) 



(D.21) 



Ho ± - r 
" 2 



M 



± 



± 



4c 



a?" f r 

± - + 1 
^2cV-^ 2 



b\r 
±—{l + 1 
2c ^ " 



Ho + 



8c2 



(D.22) 



The above equations generalize the ones in |10l HI] for general squashing factors a, b and c. 
One can check that the results of [lOl E] are recovered when a? = b'^ = Q and c = —3/2. 
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D.l The ABJM harmonics 



Let us particularize our results for the metric of the T^'^ space. In order to find the rescaling 
factors a, b and c in terms of the filling fraction u, we compare the eigenvalues of Aq obtained 
above with those obtained directly from the differential operator. It is easy to see that a, b 
and c should be taken as: 

a' = , b' = - , = / . (D.23) 

1 — u V 4z/(l — v) 

Indeed, for these values + 6^ = 4c^ and the eigenvalues do not depend on the R-charge 
r. Actually, in this case is given by: 

= + i(i±i) . (D.24) 

In order to compare this result with the eigenvalues obtained from the differential operator, 
we should relate j and / with l\ and l2- This relation is just j = /i/2, / = ^2/2 and, therefore, 
i^o becomes: 

4 1 — z/ 4 V 

which, indeed, corresponds to the values found in appendix O by diagonalizing the Laplacian 
operator. 

Let us now turn to the diagonalization of the matrix M of the Hodge-de-Rham operator. 
The eigenvalues A are the roots of the polynomial equation: 

p(A) = det (Al - A /) = , (D.26) 

where / is the 5x5 unit matrix. Notice that p(A) is a polynomial of degree five and, therefore, 
this equation has five roots. We are interested in finding the eigenvalues corresponding to 
modes that satisfy the transversality condition D"" Y(a) = 0. Since we are not imposing this 
condition in our diagonalization problem, we are also including longitudinal harmonics of 
the form VaY. The eigenvalue A for this longitudinal harmonic should be the same as the 
eigenvalue of Aq acting on a scalar harmonic, i.e. X = Hq. It is easy to determine when 
A = i^o is a solution of the eigenvalue equation. Indeed, one can check that: 

It follows that only for r = 0orz/ = l/2 (half- filling) is A = i^o a solution of the secular 
equation p(A) = 0. In these cases it is easy to identify the transverse eigenmodes as those 
orthogonal to the longitudinal one. For this reason we will restrict our analysis to these two 
cases. 
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D.1.1 Arbitreiry filling 



In this case we have to take r = 0. The eigenvalue X = Hq appears three times. Assuming 
that one of these modes is the longitudinal one, two transverse modes with X — Hq remain. 
Therefore, we get two series of eigenvalues: 

Xi — Hq , 

with Ai being doubly degenerate. 
D.1.2 Half filling 

In this case — ^ and r can be arbitrary and we get two branches of eigenvalues: 

T(±) _ ll{h + 2) , ^2(^2 + 2) 

- 2 + 2 

A?) = 1 + /i + M^i±^+M^. (D.29) 

2 2 V 2 2 

Notice that the eigenvalues A2^'' are independent of r and they are identical to the X'^^ for 
the particular case — ^■ 



E The Heun equation 



The Heun equation is the second-order differential equation given by: 

■ ■ ^ +^)^+ (E.l) 



(i^;^ V2; z — I z — d / dz z{z — l){z — d) 
with the condition: 

7 + (5 + e = a; + /3 + l. (E.2) 

The Heun equation is the standard canonical form of a Fuchsian equation with four singu- 
larities (located at 2; = 0, 1, d, 00). Recall that the Fuchsian equation with three singularities 
is the hypcrgcometric differential equation. The exponents of the singularities (which are 
the roots of the indicial equations) at 2; = 0, 1, d, 00 are, respectively, (0, 1 — 7), (0, 1 — 5), 
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(0, 1 — e) and (a, (3). The parameter q in (lE.ip is called the accessory parameter. The solu- 
tion of ( IE. II) which is analytic at z = is called the local Heun function and is denoted by 
Hl{d, q; a, (3, 7, 6; z) (see for a review). 

In some cases, the solutions of the Heun equation can be related to other special functions. 
So, when 7 = 5 = e = 1/2 the Heun equation reduces to the Lame equation. Moreover, 
when (i = 2, e = 7 and q = a(3, the local Heun function is related to the hypergeometric 
function by means of the equation: 

i7/(2,a/3;a,/3,7,a + /3-27 + l;2;) = ^; 7; 2;(2 - 2;) ) . (E.3) 

This equivalence with hypergeometric functions seems to apply in our fluctuation equations 
with E = Q. 

By performing a change of variables the Heun equation can be converted into a Schrodinger 
equation for the so-called BCi elliptic Inozemtsev model, which is a quantum integrable 
model (see, for example refs. [13 | H1 | B5]). Let us review this mapping, which involves elliptic 
functions. With this purpose, let us briefly recall some basic facts about these functions. We 
will use the conventions of [46j, where more details can be found. 

Let Ui and UJ2 be nonzero complex numbers such that Im {lo^/ui) > 0. The lattice L with 
generators 2ijJi and 2u3 is defined as: 

h = {w e C/w = 2nuji + 2171003 , n,m e Z} . (E.4) 

The lattice parameter r is defined as: 

r = — , Imr > . (E.5) 

For a given lattice L we define the Weierstrass elliptic function p{x) = p(x|L) by means of 
the formula: 

t«eL\{o} ^ ^ 

We will also define U2 as 

UJ2 = -idi - U3 , (E.7) 

and the numbers ei, 62 and 63 as: 

= p{u,) , {t = 1,2,3), (E.8) 

which satisfy the condition 61-1-62 + 63 = 0. In order to relate these quantities to the Jacobian 
elliptic functions, let us define the modulus k and the complementary modulus k' as: 

k' = , k'^ = . (E.9) 

61 - 63 61 - 63 
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From this definition one can check that k and k' satisfy the relation k"^ + k'"^ = 1. These 
moduh can be related to the lattice parameter r by means of the relation: 

K(k') 

r = .^. (E.10) 

where K{k) is the complete elliptic integral of the first kind, defined as: 

K{k) = / , . (E.ll) 

In terms of this function and the moduli k and k', the numbers Ci, 62 and 63 are given by: 

e, = ^ [Kik)r , e. = [Kik)Y , 63 = [Kik)]' . 

(E.12) 

The Weierstrass function p{x) satisfies the differential equation: 

p\xf = Apixf - g2p{x) - gs , (E.13) 
where g2 and g^ are the so-called lattice invariants, which are related to the Cj as follows: 

g2 = 2{el + el + el) , g^ = Ae^e^e^^ . (E.14) 

To convert the Heun equation into the Schrodinger equation of the Inozemtsev model, let 
us consider elliptic functions with complementary modulus k' given by: 

= 1 , (E.15) 
a 

where d is parameter of the equation (IE. II) . We change variables from the z of eq. ( IE. II) to 
a new variable x, related to z by means of the following relation: 

. = ^t^^ . (E.16) 
62 - ei 

In terms of x, the Heun equation becomes a Schrodinger equation of the type: 

H^{x) = £^{x) , (E.17) 
where the Hamiltonian H is of the form: 

H = -— + J2hik + l) p{x + CO,) . (E.18) 

i=0 
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The coupling constants k in (lE.lSp are related to the parameters of the Heun equation by 
means of the relations: 



1 



lo = (3 - a - - , 



1 
2 



= -7 + ^ , h = -6 + ^ , h = -e + - , (E.19) 
while the new function in (]E.17p and the old function x{^) (lE.ip are related as: 

^(x) = z-'-^ {z - 1)-'-^ {z - d)--^ x{z) ■ (E.20) 
The "energy" S of the equivalent Schrodinger problem is just: 

4g - A , (E.21) 



£ 



ei - 62 



where A is given by: 

^ = ^(rf+l)[(7-l)'-^(«-/3f 



2-d 



(e_l)2_(5_l)2 +27[e-c/5] . 



(E.22) 



The Heun equation with c? = 2 is a special case. Indeed, it follows from ( IE. 151) that the 
moduli of the elliptic function are 



k = k' 



1 

7! 



(E.23) 



Equivalently, one gets from flE.lOp that the lattice parameter is r = i for this case (this is 
the so-called lemniscatic lattice). Moreover, see (IE.12p . the Cj's for this case are: 



ei = -63 



2ul 



327ru;f 



2 ' 



62 



. 



(E.24) 



Let us write more explicitly the change of variables for this d = 2 case. For simplicity we 
will choose from now on the parameter Ui = 1/2. From ( 1E.16P we can write ^ as a function 
of X in the form: 



1 - 



p{x) 



in 



(E.25) 



[ni 

The function p{x) in ( 1E.25P is a doubly periodic function in C with periods {2ijJi, 2U2) = (1, i). 
When X is on the real line p{x) is real, it has a minimum at x = 1/2 (p(l/2) = ei) and 
p{x) — )■ +00 at X = 0, 1 . This means that z runs twice over the interval [—00, 0] when 
X G [0,1]. Notice that, to parametrize the hanging brane configuration of sections [6] and 
[HI the coordinate p defined in (18. 2p must run twice over the interval [1, 00], which precisely 
corresponds to taking x G [0, 1]. 
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It is also interesting to write the inverse of the change of variables (1E.16P for the d = 2 
case. By differentiating ( IE. 161) and making use of ( IE. 131) for (72 = 4e^ and = 0, one gets: 



^ = ±2,/^i^z{l-z)iz-2) 



(E.26) 



where the + (— ) sign must be taken when x G [0, 1/2] (x G [1/2, 1]). This equation can be 
integrated in the form: 



X 



1 ± 



2V27r 



dt 



t{l-t){t-2) 



(E.27) 



where we have already taken into account that x = l/2for2; = 0. 
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